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ABSTRACT 

The cooling load for a railway coach has been computed 
from the governing ec^uations formulated on the basis of heat 
and mass transfer analyses, A generalized computer program 
has been formulated for the same for three main routes of 
Indian Railways for which evaporative cooling can be used,' 

The results have been obtained for the coach movins with 
a given velocity for various places with the known latitudes, 
longitudes and outside summer design' conditions. 

The non-regenerative and r egenerative evaporative 
cooling systems have beon considered. They have been 
compared with each other in terms of tonnage and the comfort 
conditions within the coach. The latter is found to be a 
better choice from comfort point of view,' 

The higher inside condition without violating the 
comfort requirement is desirable as it can be achieved from 
the evaporative cooling with air flow rates ranging from 
10,000 to 15,000 kgAi.' Fbr required comfort conditions, the 
computed results reveal that even 5 -ton air conditioner 
will serve the purpose as compared to the existing lO-ton 
air conditioners used in railway coaches,' 



CHAPTER 1 


1.1 INTRODUCTIOH 

Ihraporatj-ve cooling is the oldest known application 
of air-conditioning principles. In the pre scientific era, 
it flourished as a natural process of cooling mainly in the 
hot-dry and lightly populated desert areas. Even in the 
scientific era, it did not get much attention as is evident 
from the fact that as late as 1958, the official American 
Society of Heating and Air-jconditioning Engineers* Guide 
comprising 1272 pages? devoted only half a page and that 
too without tables or data. But, it started gaining due 
attention in the modem scientific era, only after the 1960*s 
when the cost of energy started escalating due to severe 
energy crisis. Hence the mechanically 'refrigerated 
conventional air-conditioners became very expensive beyond 
the means of average income -group people. Thus the former 
was looked to be a promising cheap alternative method of 
air-conditioning. This is supported by the fact that the 
ASHRAE Handbook [ 1 ] devotes a separate chapter on 
evaporative air cooling and related equipment. Howef-er, it 
did not receive much scientific impetus then as it was 
considered to be feasible only in particular areas and is 
impractical for other areas ,like the coastal regions. 
Unfortunately^ it was not realised that the sea itself acts 
as a natural evaporative cooling device as is evident frcan 
the temperature of the mainland compared to only about 35°G 
of the sea coast. Hence, ,ln the coastal areas , the human 
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comfort conditions can be maintained by high rates of air 
circulation over the body. 

Prehistoric evaporative cooling started around 
3000 B.C., when in India, evaporation was even used to make 
ice [ 2 In rural India, the people still use the 
ancient cocos tatti of great effectiveness. In olden days 
doors were replaced in summer by tatties, frameworks similar 
to screen doors covered with dried khuse -khuss grass which 
were kept wet by coolies, Ihe most popular is the Indian 
evaporative air cooler, known for over $0 years where 
instead of door^a wheel like framework covered with khusS~ 
khuss grass is revolved and kept moist by coolies, Liko 
tatties, the grass lends a pleasant odour • Another classic 
example is the use of earthen pots for cooled drinking water. 

These days evaporative cooling is gaining much 
importance!, Several workers have tried to get optimum 
performances of desert coolers [ 1 based on pad 

density, w^ter spray over the pads etc.' The evaporative 
cooling has also been used in cars and in buses [ 6 1 because 
of cheaper system. A proposal is under way for uae of 
evaporative cooling in Indian Railways [7 1* 

1.2 THSORI of EVAPOEATITHl COOLING: 

There are two types of evaporative cooling ftiirect 
and indirectt, The principle underlying the former is 
conversion of sensible heat of air into latent heat of 

vaporization of water added directly into the ‘air* In the 



latter, the air is cooled by the evaporation of water not 
contacting it, so that the sensible cooling of air occurs 
without increase in its moisture content [83* 

The exchange of sensible heat for latent heat under 
the isolated conditions results in the lowered, temperature 
of air and the air becomes eventually saturated at the 
water temperature. This process is called ’adiabatic 
saturation’ since there is no external heat exchange between 
the ;^stem and surroundings. The practical applications of 
this principle is found in pump-equipped drip ~type 
coolers, slinger and rotary pad-coolers, most commercial 
air washers having spray -chamber or capillary design and 
some textile -mill cooling systems. 

In practice however, the water usually gains some 
external sensible heat in the sump tank, pump and piping. 
The make-up water entering the sump to replace the evapora- 
ted and lost water adds heat. Other sources of heat 
addition include circulatory friction, heat transfer from 
the surroundings and possible solar radiation. Thus most 
adiabatic saturation’ in evaporative cooling is merely 
a close approximation, with considerable evaporation 
occurring to cool or recool water. When the circulating 
water is considerably warmer than the air wet -bulb 
temperature during initial contact, the process resembles 
that occurring in most cooling towers for cooling warm 
condenser water. Air and water are jointly cooled and not 
the air alone. This occurs in direct evaporative cooling 
systems having ’ once -through' or purapiess use of water, 
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as in many small drip-type coolers, some spray -chamber and 
capillary air-washers. More water is consumed and the 
cooled air is both waimer and more humid than in true 
adiabatic saturation,' 

There are limits to the cooling achieved by adiabatic 
saturation. The amount of sensible heat removed equals 
the latent heat required to saturate the air with water 
vapour* The cooling possibilities of adiabatic saturation 
vary inversely with the existing degree of humidity of air,' 

Most commercial direct evaporative coolers deliver air 
cooled about 70 to 95^ tw-fards saturation. This is only 
50-70 percent for inexpensive drip-type coolers, but for 
slinger and air-washer types, this amounts to over 90 
percent of perfect humidification, A ratio known as 
saturating or saturation efficiency (n ) is used for 
rating the performance of such devices. It is given bys 

^ = - fab I . C1.1) 

- ’^wb 

where = entering air dry -bulb temperature (°C),' 

= entering air «et bulb temperature (“O.' 
and T =* the dry-bulb temperature of air leaving 
the system 

However, the performance of evaporative coolers 
will be badly affected if heat enters the process from any 
source other than the air being cooled,' In view of this 
it is recommended that evaporative coolers be (1) Ipiqated 
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in the shade whereTsr possible, Cii) have the coolest 
possible water supply and (iii) receive the coolest and 
driest air,' 

l,-3 USB OP Wi-?0RATIVP1 GOOLIFG IH RAILWAX COACHES 

The Indian Railways also desires to use evaporative 
cooling , especially in the three -tier sleeper coaches because 
the passengers feel very uncomfortable during summer 
season. There is considerable scope of evaporative cooling 
for trains plying in the areas aiway from the coastal regions. 

Ife met the authorities at RDSO 'Uicltenow In JUly 1981if* 

They supplied us various drawings and relevant data* W© 
also surveyed the existing air-conditioning systems of the 
railway coaches.' As such, they do not possess any standard 
method or data for cooling load calculations [93r’- They use 
lO-ton conventional air-conditioning unit divided into two 
separate units of 5 ton each. The cooling of air is done 
from both ends. The air is distributed uniformly through 
a central duct on both side/^* 

So we carried out the preliminary analyses and the 
results published recently [ 10 ], Further, calculations 
were carried out on the use of evaporative cooling for 
railway coaches. The comfort conditions were estimated 
using Fanger* s comfort eq.uation for temperatures below 
28 It is found that the same lies in the comfort zone.' 
Use of evaporative cooling with regeneration reveals that 
it can give the desired comfort conditions at the cost of 
additional equipment. 
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1,'4 PHSSBNT GOOTRIBtJTION 

A computer program has been developed to calculate 
the cooling load, temperature distribution Inside the coach 
and other parameters of interest at any place for a given 
time and velocity of moving coach with known outside summer 
design conditions. Also, the effect of the mod© of inside 
convective heat transfer on the overall heat transfer rate 
has been studied. The heat and mass transfer analyses 
have been carried out for two separate cooling systems for 
evaporative cooling in the railway coach and their usefulness 
ccHnpared in terms of energy saving and achieved comfort 
conditions inside the coach. The effect of various coach 
parameters on the tonnage and temperature distribution 
inside the coach has been thoroughly investigated. 

It has been found that the regenerative system 
yields better comfort conditions but involves more number 
of components* The computed cooling load is found to be 
barely half of the capacity of existing air conditioners for 
the railway coaches. 

It has also been found that the water requirement for 
evaporative cooling for 6 -hour servloe duration is found 
to quite reasonable for being carried by the train," 



CHAPTER 2 

MATHEMATICAL FORMHLATIOH 

2^1 DESCRIPTIOH OF THE SISTEM 

The air-conditioning of the railway coach is done using 
mechanical refrigeration system. The compressor and condenser 
are mounted underneath the coach. The air-cooled condenser 
having a blower is used to cool compressed vapour. The 
condensate is supplied to the cooling coil mounted at the 
ends of the air-con (31 ticned coach. The air is distributed 
through a constant area duct to both sides of the coach. 

In the present analysis, the evaporative cooling 
is envisaged to get t^perature in the range of 28 to 34 
The air is supplied from one end and is exhausted to the 
other end. The evaporative cooling units are mounted in the 
space of the vapour-compression machine. Even water-tank is 
accommodated in the same space. This is done in order 
to avoid further complexity,' 

2.2 MATHEmTICAL TREAT MEHT 

Figure 2,1 shows a 3 Tier broad gauge CSC 1640 
II class sleeper coach with dimensions, A differential 
segment dx has been marked showing heat and mass transfer. 
Then, for the control volume we write the mass and energy 
conservation equation as: 
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GOVSRJJIHG IQUAIIONS 
2,2*1 Mass transfer: 

The conservation of moisture component, is: 

- tt t - ^ 

dw 5= w , A^. dx <2*1) 

where A- = cross-sectional area of the coach (m^l, 

= mass flow rate of air ( kg/h ).‘ 

Ht 

w = moisture release by persons per unit time per 
unit volume of the coach (kgw/h-m^), 

w"’ = np/(A^ X ) (2.2) 

where = moisture release per person per unit time 
( kgw/h-per son ), 

Up = no, of passengers in the coach* 

= length of the coach (m). 

On integration of Eg., (2*1) with boundary condition 
(x=0,w = W2) one gets ; 

w =5 (w**A^/iftg^ ) X + W2 (2,3) 

where W2 == specific humidity of air at inlet to the coach 
(kgw/kg dry air), 

w = specific humidity of air (kgw/kg dry air)* 

2,2,2 Heat transfer: 

The heat transfer equation for the control volume is 
written as: 

dh = 


( 2 . 4 ') 
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where = occupant heat release inside the control 
volume (kJAi ). 

dQg = heat input to control volume c3x due to solar 
radiation (kJA). 

(IQq = heat input to control volume dx due to 

conduction through walls, floor and ceiling 

(ki/h ), 

= total heat input through structure due to 
conduction and solar radiation to control 
volume ( kJ Ai) • 

The specific enthalpy of moist air is expressed ast 

h = c^^ T 4 - w ( 2?01.4 + 1.884 T ) (2.5) 

pa 

where T = temperature of air inside the coach (°C), 

Using Eqs. (2.3) ahd (2,'5), we have 

h = Cpg_ T + [(w"’ ) X + w^ 1(2501.4 + 1.884 T) 

( 2 . 6 ) 

and using c„^ = I.004 we have 

g = { 1.004 + 1,884 1*884 w"\/\) X. + 

(1.884 w"'Ao/aa 2501.4 w"'a^/iI1 ^ (2.7) 

The heat generation is written as: 

dQj^ = dx (2.8) 

lit 

where q is, the occupant heat generation per unit volume 
of the coach (kj/h-m^). 


(2.9) 
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where H = occupant heat release per unit time (lcJ'/h~person). 

Jr 

The structure load dQgg, can be determined after 
making the following preliminary calculations: 

(i) 4. ( Area of surface, m )i 

A. = f..dx (2,10) 

J J 

For side wall (j =1); 

A-i [( n^ ^^“^^1 ^ dx = f^.dx (2.11) 

where, = height of the coach 

j 

n^ = no, of windows on one side of the coach, 

f 

= Area of one window (m'^), 
ceiling ( j =2 ) : 

Ag = [(P^- 21^- L^)/23dx = fg.dx (2.J2) 

where P = perimeter of the cross-section of the coach <m), 

= width of Ihe coach (m) ' (Fig, 2»2), 


For glass window ( 3 = 3)5 


A^ ^ C ( ^w^ ^ ~ ^3* 

( 2 . 13 ) 

For floor ( 3 = 4)1 


= (L^/2).dx = f dx 

(2,14) 


<ii) (%)gj (Intensity of total impinging solar radiation 
on the surface, kJ/h-m^); 

(a) h (Hour angle, degrees); 

Local Apparent Time (L,A,T, ) of the given place [ll] 
is related to the local Mean Time (I%M,T,) as ; 
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L. A.T, = L, M,T, 4- Equation of time vdth given by 

L, M*T. = I. S, T. +• correction where I, S*T, - 3hdiaii S'tandard 


Time 


and correction = k - (L^) gjjj](minutes) 

with ” Longitude of the place in degrees# 

(Lo)gm“ Longitude of standard meridian in degrees 


m -^s 


For May 21, equation of time = +3 37 


L.A.T. = UlUT. + 3“ 37® (2#15) 

The hour angle (h) with L. A*T. in hours is given by: 

h = 1? (L. A.T - 12.00) (2#16) 

(b) (Intensity of solar radiation normal to sun^s rays 
at the location, 

The solar altitude angle ( ^' )[ ip] latitude of the place 
(1) , declination angle (d) and the hour angle (h) all in 
degrees., are related as [12]i 

sin = cos(l) cos(d) cos(h) + sin(l) sin(d) (2.17) 

We have 

Ipjj = A/exp (B/sin^) (2,18) 

where A and B are constants depending on the time of the 
year [I3], 

(c) (Intensity of direct solar radiation, : 

from Fig. 2*2, one gets: 

cos # = (sin|p sin 1 - sin d ) / (cos^ cos 1) (2*19) 

where § ~ solar azimuth angle measured from south 
For solar noon (h = 0), f = 0. 
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Pc s p«rirTiftt«r of cross section of the coach (m) 
Af *Area of cross-section (m^) 


inNErt WALL 
SURFACE 


Fig. 2.2. Cross-sectional view of the inside wall of 
railway coach. 


SUN 




Fig.2.1 Var!Ou:v solar angles. 



( 2 . 20 ) 


For sidewalls and window ( ;j = 1 and j = 3)j 
cos CSgj^ = cos P 
where e . = angle of incidence of solar radiation* 

S J 

/ 

Vo-i ~ surface solar azimuth angle measured from south. 

If = surface azimuth angle measured from south, can 

he obtained from the known values of # and for various 

s 

train routes (Table 2,1 ): 

TABI^ 2,1 

VALUES OF W AMD V FOR DIFFERENT TRAIN ROUTES <Fig, 


Train 


Orientation East facing side 


West facing side 


Route 

of Cbach 

of 

coach 


of coach 


’*1 

^“(a. 


SouthJSast 

NW-SE 

13?° 

»-»s 

3 

45° 

s 


South 

N-S 

90 

mm 

S 

§+T 

S 

90° 

S 


South-West 

NE-SW 


m 

1 j 

s 

135° 


If / 

sj > 90° i. e. cos 

< 0, the 

surface is in 



the shade and the intensity of direct solar radiation 
impinging on the surface ( Ij^) = 0, 

For celing (;} = 2 ): 

c"os (Ogj) ~ sin^ (2.21) 

Here the ceiling has been assumed to be approximately 
horizontal in spite of its curvature: 

Thus, , 
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(d) Intensity of diffuse radiation (kJAi-a^) 

For sidewalls or glass window (;) =1 or ^ =:3)s 

where C = diffuse radiation factor. 

= sky diffuse on vertical surface 
sky diffuse on horizsontal surface 

where is given by 

Ig^ =0.5? + 0,437 cosCOg^) + 0.313 cos^COg^) for 

cos e„. > -. 0.2 (2.23) 

S3 

For cos 0g^ ^ - 0.2 Yg. = 0,45 [H]. 

For ceiling ( 3=2); 


^^d^sj “ ^ ^^ss^sj 


<2. 24) 


where (F^ )„.= angle factor between the surface and the sky. 
ss S3 

If ( ^)_4 - tilt angle of the surface measured from the 
s J 

hori25ontal, then; 

*^*■85^53 = < 1 + OOS )/a = 1.0 


• • 


*83 = ° 5 


(e) (It,).,^ (Intensity of reflected radiation, kJ/4i,4!i^) : 

^ S j 


where 


Ii.TT P ,-*» (F„^), 


gr 


tH ^^gr ''"sg^sj (2.25) 

= total radiation ( including sky diffuse and 
direct) falling on the ground (kj/b-jn^). 

= reflectivity of ground. 


(P ) .= angle factor between the surface and the ground. 
Sg ’ S3 

We have .. , 

^tH ~ ( C + sin^ ) 


< 2 . 26 ) 
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= [ 1 ~ cos 6g. ]/2 


(2.27) 


For sidewalls and window- ( j = 1 and j ~ 3) 


^ . = 90 ° 
sj ^ » 


» 

K ♦ 


(F ) . = 0.5 

sg^S3 


For ceiling ( j = 2) : 


SD 


0. 


f 

« * 


(F . 
sg S 3 


0, 


Hence no, reflected radiation falls on the ceiling. 


From the above the total radiation (l^,) impinging 
on a surface is found as: 

For sidewalls and ceding: 

= Vs3 

For floor: 

(Ij) = 0 since negligible solar radiation is assumed 
to impinge on the bottom of the coach, 
(iii) h^ ( outside convective heat transfer coefficient, 
kj/m^ - h - °C): 

Assuming a transitional Reynolds number of 3' x 10^ 
and a train velocity V’C km/h) greater than or equal to 
60 kmAi we have 

^cr ^1 (2,29) 

so the whole flow can be assumed to be turbulent. Hence 

= 0.0292 <Re^)°*^ <Pr)^/3 Cz*30) 

for forced convection turbulent flow over a flat plate [15] 
Using Iq,' ( 2 * 30 ) and taking the properties of air at the 
mean coach outside temperature lof 38 °C , we get: 
h^ = 0.02lf3 (V X 1000)^*® 


(2.31) 



18 


Xjjj = mean distance (m) of the control volume from x = Or 
The above eq# (2*31) holds good when the supply of air 
is in the same direction as the velocity of the train* If 
they are in the opposite directions, we have 

hjj = 0,02l<-3 (V X 1000)^*^ (L^ - (2*32) 

(iv) R . ■ (Besistance of structural component, ra^- h - ^C/lcJ-) 

J 

The values of resistances, of various structural compo- 
nents of the coach viz* sidewall, glass windcw , celling and 
floor have been calculated according to their individual 
structural details given by the drawing [i^] and have been 
found from the air gap resistances [17] ass- 

Sidewall j = 0.0560 - h - ^C/kJ. 

Ceiling : R^ = 0,0704. m^ - h - ^C/kJ, 

Glass : R3 = 0V00i27 - h - °C/kJ which is quite 

Window negligible. 

Floor : \ = 0. 1139 m^ - h - ^C/kJ^ 

(v) (Inside convective heat transfer coefficient, 

kJA-m^) : 

Two cases have been considered to study the effect of 
internal mode of convective heat transfer on the net cooling 
load • They are: (a) CJiSE I: Forced convection for celling 
and free convection on other sides of the coach and Cb) 

CASE II t Forced conyection on kll sides of the coach* 

(a) CASE It Forced convection for ceiling and free convection 
on other sides Of the coach* 
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For ceiling ( j = 2) s 

Assuming a transitional Reynolds number of 3 x 10 
Xqj, is calculated frcan the forced convection velocity of 
air flow under the celling . If JSj.. is the net flow 

area (m ) then 

%c “ \ ^ 3600) (m/s) (2.33) 

where v = specific volume of moist air C m^/icg of dry air) 
whose value is taken to he 0*87 m^/kg dry air# 

Olhe values of estimated frcm Bq# (2.33) ior 
different air flow rates range from 0,5 - 2*0 m/s and the 
corresponding x„_ values have been found to b e greater 
than half the length of the coach Hence the flow 

beneath the celling is partly turbulent and jartly laminar 
and an average film coefficient (^^av^sj ^ found 

from the following [ig ]4 

Nuj^ = 0,036 (Pr)^/3 [(Re^)^*^ x8700] (2.3l+> 

from which we have 

0.'0302 x 3600)°*®(I,j_)°‘®. 57.6]A-i 

-“O C2.35) 

Other sides ( j / 2 ): 

Since free convection has been assumed, Grashof 
number is calculated for all the other sides# 

For airs 

(Gr*Rr)g^ = 1,6 X 10“^ 1^^^ ( 

where 1^. - characteristic length (m)} # 


( 2 . 36 ) 
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= teiaperature difference between that of the 

S J 

surface and the flowing air C®C), 

For sidewall ( j = 1 ) : 

Ig^ = and 10^ < < Gr*Pr)g^ < 10^^ (2*37) 

From Eq, (2.37) it can be deduced that free convection 
flow over the sidewall is turbulent* Hence, 

(h )g^ = l}..74 ) ^ ^ vertical plates) 

( 2 * 38 ) 

where 

(hi)gj = inside heat transfer coefficient (kJ/h-jn^- °C) 

=5 temp* difference between inside wall 

surface and Inside mean air temperature (°C) 

( ^^vi^s^ given by: 

^ ^\i^s3 " ^^wi^sj " (2.39) 

where 

(T^)g^ = inside surface temperature of the coach (°C), 

= mean inside air temperature of the coach C®C), 

Since the coach has been divided into 23 elements of equal 
length, i varied from 1 to 23. Hence 

^Vi = ( >/2 (®C) (2*40) 

where = Temperature of inside air at the i^^ section (^07.: 

^i+1 ~ temperature of inside air at the (i + 1)^^ 

section (°C), For i = 1,2,. ..23 

Fbr glass window ( j » 3) f 

igj ~ ^gW ^ ^ 


(2*41) 
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where h = height of the glass window Cm), 

5 V 

Hence free convection is laminar and 

(2.42> 

( for vertical plates): 

For a known = 0,600 m Eq, (2. 1+2) turns out to he 

(hi)^j =5.80 ( (2.1,3) 

For floor ( ^ ~ k) ’• 

^sj = and 10^ < (Gr.Pr)^^ < 10^^ (2,1+4) 

Hence free convection is turbulent. However the floor 
is equivalent to a heated plate facing upward* So the 
appropriate equation for this is 

(hi)^j = 5.47 ( <2.45) 

For j 7 ^ 2 , ( determined using ‘Newton-Raphson 

method*. We have: 

C2.46) 

where U . = overall heat transfer coefficient (kF/^-m^- ^C)'i 

Sj 

(Tgoi^gj = sol-air temperature (®C) the evaluation of 

which is explained In Cvl),‘ 

T = temperature of air inside the coach (®C)* 

Also, == l/l(l/h^)+ + l/(h^)gP (2,47) 

Using Bq. (2*46), we haves 

«80l>sj - * = < V.l < ^ <2.48) 

Putting ~ a»<i using Bq, (2*47) |W« have: 
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^'^sol^sj “ ^ "^wi^so 


\i^sj ’’' ^ ^■wi^sj ^^3 


•30^30“ - - ^ S‘ Wsj' ■ ^ 


^ Wr- " 

Putting - T = u’ and C^(R^ + lAi^) u”, w# have (2.49) 

(2.50) 

J 

Eq. C2,'4'3)is of the form: 


n m.+l 

^ ^ ^ ^ -u'=0 


f('^Wsj=° ^2-5« 

and using Newton-Raphson iteratioe technique to determine 




"’^vP.3 


where [f ' ( ^ " T~ — C 

^ d(A T , ) 


wi^sj 


f' ^ ^ ^ ^ ’^wPsS 


(2.52) 
and is given by 
ia3 


(2.53 


The ahoi^ iteration is of second order. 


(b) CifiB II: Forced convection on all sides of the coach: 

For ceiling; 

The same eq* (2.35) used for case I is applicable. Hence 
(hlavhs “ ^-^fo * 3600)°*® (Lj^)°*'® _ 57.6]^^ 

(2.54) 

For sidewall and glass window: 

^cr ^ ^ ^ transitional Reynolds number of 

3 j 00,000 wMch means that the flow is laminar throughout the 
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surface, 

, Nuj^ = 0,664 (Pr)^^^ C2.'55) 

or, (h^aTPg^ = 0,237 

For floor: The flow is partly laminur and partly turbulent 
as in the case of ceiling and again we have 

t 0-0302 CVj^ X 3600)°* 3 57. 6] /L^ 

(2.56) 

^ ’^sol^s3 = - 

For ceiling and sidewall * 


<Isof S3 = "db + “^3 ^ Vs3/ ‘'o 


(2.57) 


where, 


solar absorptivity of the surface. 


ForPioors *; (It)sj= 0! (IgofsS^’^db 

For glass window: 


(a?e: 


(T 


sol^ s j 


where 





ds 


rs 


~ ^^s’^Fs '^6.3 ^^d^sj 

+ %s (2,59) 

= solar absorptivity of glass for direct radiation. 

= solar absorptivity of glass for diffuse rac’lation 

=5 solar absorptivity of glass for reflected 
radiation. 


F_ = shading coefficient assumed to be 1,0 due to 
s 

negligible shading over window glass, 

Ds» *^“rs determined from the following: 
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Assuming coratsion window glass tjf thitjkness of doulDle- 


strength c[uallty, we have: 
-K t. 

a_ = e 


a 

G 


( 2 . 60 ) 


where ~ absorption coefficient of glass* 

t: 

K - extinction coefficient [ fq]. 


tg/(/l -"sTn^ ^ 


(2.61) 


where n = index of refraction for glass* 

& 

= incident angle for radiation component ( ° ). 

w S 

Direct radiation component; e„„ = e„. (2,62) 

cs S J 

D ( c = 1 ) 

Diffuse radiation corapcaaent ; e^g = 60°, A mean angle of 
(c = 2) 

incidence of diffuse radiation has been assumed [20 3 
Eeflected radiation component? 6^^ ~ 50** assumed [g’l ]. 


Now T 


(1 ~ r )‘^a 
^ cs g 


cs 


1 - r 


cs 


a 2 

§ 


(2,63 ) 


where, 

T ^ = transmissivity of glass for single radiation 
c s 

component,' 


r 


cs 




cs 


specular reflectivity for single radiation 
ccanponent, 

^ ** ^cs - ^ ^cs^^ ■“ ^cs ^ (2*64) 


and 



where 


“ absorptivity of gXas-s for single radiation 
component.' 

Using as reflectivity of glass for single ra(3iation 

component, one gets pQg = l - ^ as 

P + + “ I,' 

cs cs cs 

Direct radiation component: ^ Ds ~ *^ls 

Diffuse radiation component: -r^g = T2 s ’ ds ~ ^2s 
Reflected radiation component rs ~ ^3s 

To calculate the amount of radiation ( 0 .(kJ/m^-h)) 

0 

transmitted through glass window we user 

^s “ ^^s "^Ds \s ^ds ■** *^rs ^rs^ (k:JA-iii^)<2.65) 

F = shading coefficient assumed to be 1*'0 due to 

s 

negligible shading of window glass* 


Finally, we haver 

- Z £ U . [(T . T] A, + 

3^1 s=i ^ sorsj j 

3 2 

"SI <^0* (2.66) 

3=3 3=1 ® 3 

Using Eqs* (2.4), (2.8), (2,9), (2.66) , we have: 

4 2 

!n^ dh = <1". ^ X 


4. 


0 .f ..dx 

3 5 


(2*67) 
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From Eq, (2*67) we have: 


4 2 



2: 

s=l 




3 

2 


4- 

X 

X 



0=3 

s=l 


(2.68) 


Also using eq« (2*7), one getsr 


fii. 


dh 


dl 


"a dx ~ ^ 1.884w^) g +(1.884 v'\)x.^^ 

+ (1.884 w”t A^)T + 2501. 4w’‘’^ (2.69> 
Equations (2.68) and (2.69) yield: 

(1.004 + 1.884 w^) + (1.'884 X ^ 

+ (l,-884 w»» A^)T 

H t H t ^ 

+ 25’01.4 w -^o “ ^sj^j^^sol^ s3 


0 -f .,dx 


j=3 s=l 

4 2 


3=1 s=l 




(2. 70) 


How putting 




a = 

1,004 

+ 1,884 w^, 

it 

< 

HI 

4 

2 


4 

2 

X 

5 

^s3 

X 

X 

3=1 

to 

3=1 

s=l 



3 

2 




4. 


a 



II 

*nj 

s=l 



U I 


°s3 


Bq. (2,70) reduces to the form; 



+ bx 


(2.71) 


a fn 

a 


42 

dx 


41 

dx 


+ bT + 2501, 4 T -h s. 


.*, dx " 2501,4- e (2,72) 

Further putting, = (b + m ) and 5*2 ” ® -2501,4 e) 

¥e lia\7e, 


( afc^ +ta ) g +Fj_ 1=^2 

(2.73) 

r ♦ e • , 

(a &^ + hx) = F2 - F^ T 

< 2.74) 

1. e, , 

r dx ; _ I dT 

J a fQ^+ bx''- J - Fi T 

(2.75) 

i. e. 

log^( a fcg^ + bx ) = 

(2,76) 

where 

C = constant of integration 



Using boundary conditions x = 0, T = T^, where T 2 = temperature 
at inlet to the coach (*^0), we have 

C = -|log U loge (*'2 - ^^ 2 ^ <2-77) 

Using Sqs* ( 2 . 76 ) and (2,7'"'), ve have; 

ilogj a r ta ) = _ logg(F 2 -^l I) + a \ '> 


i# '©• 


1 , , a \ 






-^1 

Tfr -m ITT' 

, 1 . / ^2 ^ *^1 “^2 \ 
lOg^ ( ^ ) 


F2 - F3, T 


(2.78) 

(2.79) 


■K^ F^/b 

• (F^ - F, T^)/(F^-B’, T) = (1 + (2.80) 

. . 2 ^2 2 1 a Bi^ 

Eq, (2.80) is the governing equation for the inside temperature 
distribution of the coach. 
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2*3 EmUATION 0? COACH INLET AND EXIT GONDITIQNSt 


(a). CASE -A ( without regeneration) i 

2.5 

Figures/(a) and (b) show the non -regenerative cooling 
system and the process during humidification in the humidi- 
fier, Assuming as evaporative cooler efficiency (r? ) ot 0,75} 
we haver 


f 1 > +’1 ^2.81) 

where (°C) and (®C) are the summer design conditions 
of a Particular place for a particular month [ 22 ] 

Using Carriers equation, 


Pv = Ps - tP - Ps‘\b^^ 


(T- - T , ) 

— ^ iHS (2,82) 

151+7 -1,1+1+ T ^ 


where the function p <T) is given by [ 23 . 3 

o 


log - ,a 2lj ^ ^ [7.21379 + (1.1520 X lO-^-i+,787 X 10*^ 

Ps(T) ^ 

(T + 273.16))(T-210,0)2)'3c a.) 

^ (T+273.16) 

( 2 , 83 ) 

where T is in °C and p„ (T) is in bars.' 

s 

For outside air the specific humidity w^(ls:gw/kg of 
dry air) at the design condition is found from : 


Wg^ = 0,622 /(p -p^) 


(2r 81+) 


~ relati.ve humidity of ambient air =? 

( 2 . 85 ) 

For inlet conditions to coach we have 


Wg = 0,62a(i^)^/[ P-(p^)j, ] 


( 2 . 86 ) 



RAILWAY COACH 


Li *21. 337m 


D8T C*C) 


Rg. 2.5. (a) Schematic representation ot non-rcgcnerative evaporative 
cooling system. 

lb) Process of humidification of - regenerative cooling. 
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4 “ ^^ 2 ^ ( 2 . 87 ) 

where (p^)2 is given by 

(P^)2 = P, (T^fe) - T„b > 

( 2 . 88 ) 

and $2 “ relative humidity of air after humidification in the 

humidifiers 

At coa,ch exit vre have 


w 


e 



A 


c 


tn 

"a 


Lf + 


<2. 89) 


Also, = 0,622 (p^)g ~ /(p-(p^)q)* 


(2.90) 


From Sq., (2.90) we have 

(p ) = JL2!i 

0.622 

6 

and 


(2.91) 


(b) (jASE_B (with regeneracion) : 

In this case, a fraction of exic air from the coach is 
recirculated through an indirect contact regenerative heat 
exchanger having effectiveness s ( Fig, 26(a)). The ambient 
air gets cooled from state 1 to state 1* without change in 
specific humidity w_ (sensible cooling process) before 

Q. 

entering the humidifier. The air is then humidified having 
process line 1* 2 ( Fig. 26 (^))» inlet air to the 

coach at state 2* leaves it at state 3, a part of which is 
recirculated.’ The above process is repeated until steady 
conditions are achieved.' 



AMBIENT 'AIR 
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■''"c-’ Vl-'j.’ ' if :> 


n iicm vij: £X!7 Ar? 

.1 f mm 1*<!I i‘ I 


I 

iHUMiOlFiES'i ■■ 
. ! f 


RAHWAY COACH 

mf 9 M \ 

‘2>'^2 ! 

l’^-- Ljx 21.337 m - 



3 


Ta, wa 


\r£JSch~Q ■ 

;';ea 7 
EXCHA^ 




r •> ': i, (a' ‘A'D/jts^n^ot'O.n of "*®c«rC'>ra!ivs evaporofivc 

LcC.AC ,;/s;err. . 

^’■ncpss c’ hu’Tiidif icQtion for f'Zgensra'iivc’ coo*.sng 
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lAfe have: 

\rhore n is ths number of iterations under the following 

Imposed const rants: 

(1) The temperature difference between the ambient air and 
the recirculated air should be greater than 3*^0 for 
useful regenerative cooling: 

’^db - > 3 °C C2.93) 


(ii) For steady state: 

l^^a^n+1 “ ^ ^ (2*94) 

The method of estimation of (T^)^ has been described as 

1 n 


under : 

= ^2.95) 

Also, (wpy = 0.622 /(p- (Py^n ^ (2.96) 


From Eqs* (2*84), (2*91) and (2*-92) we have 




(2.97) 


If is the wet bulb temperature of air (^C) after 


sensible cooling in the heat -exchanger, then using Carrier’s 
equation we have 

(ft ) = P (Tv)- Hk-JL- (2.98) 

5s Wb -'n 1547-I.lf4 (T„t,')„ 


IS o (T ') t"-V"wb'a(Ia>n-f^b'U 

i.e, Pg ^ 


151 ,^ 7 - 1 . 44 ( 1 ^^ )„ 


-(p ^ ) =0 

'^v n 


(2.99) 
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vhich can be put in tbe form 

= 0 ( 2 . 100 ) 

This is solved for '^wb ^n using Newton-Raphson method 

as detailed in Appendix b«' 

With Icno^m (!„>,*) and be found fr<an 

wD an' in 

the following: 

(Tl)j, = (Tg^)^ ( 1 . „ ) + n (T^b'^n 

For state 2 we have 


[(P ) 1 - -D (T S • 

. V2*^n ^n “ 

Fi'om eq, (2.102) we have 


1547-1.-44CT^^*)„ 


<2*10^ 




(2.103) 


and, 


(§2 


*“ t^Py^2 ^n ^^i^n 


(2.104) 


From the known above inlet conditions to the coach, and by 
using the governing equations of heat and mass transfer 
for the coach, a new value of (Tg)jj was found which was subs- 
tituted in Eq* (2.88) to get another new value of 
This iteration process was continued imtil Bq. (2.94) 
is satisfied within the preset limit* 

For coach exit putting x = In Bq,' (2.'3)| we have: 


Also, 


(w ) 


0 , 622 ^^%^ e^n 


P- [<P7>e^n 


(2.106) 
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Fj?cta E<i, <2* 106) we ii«ve« 




iiVn, 


0 , 62 t. + CWg)^ 


»e ("5 = [(P^)g]n /Pg (Vn 


(a.107) 

(2*108) 


Thereafter the following quantities were calculated as a 
final result: 

No. of iterations n. 

After iterations the final values are: 



» 



CHAPTER 3 

COMPUIAHOU OF VA.iIOUS QUAUtlTlES 

This section deals with the computational details 
used to determinet 

(a) Cooling Load (TR) 

(b) Ambient and Coach Conditions^ 

(c) Water requirement for the humidifier 

Cd) Air washer dimensions of the humidifier and 
(e) ♦ Sensible heat Conversion Rate (SHCR)* and *tJnit 

Sensible Heat Conversion Rate (USHCR) ’ 

The computational methods for the calculation of 
the above-mentioned parameters are given below: 

(a) Cooling Load (TR) : 

As the overaell heat transfer coefficient Is 
implicity related to x the Integratioi;;! of the goTeming 
equation cannot be done for the entire leng'tte^ 

Hence finite element method has been used# The coach has 
been divided into 23 sections and the total heat load 
^ ^TOTAL through each section has been calculated 

and summed up to find the grand total heat load 
QjoyAL^l’^'TAi) for the coapb. Mathematically* 

For. J = 3; 

2 

+ <?, (I.i/23) 


(3.1) 
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For j 9^ 3: 


( ~ i'-f ”* ^7T» ^ 


S=1 


si ^^-^soVsi • -"m 


(3*2> 


where = me^-n temperature inside the section which 


1.. given hy 


® ^ ^i+1 


(_^.3)' 


The heat release ^ (kJAi) , is found as: 

«R = <1'" ^ <.\/ 23 ) 


(3«4) 


The total heat load is calculated as: , ;., 

^ ^ ^TOTAL ^i ^ ^ ^ 

Henoe the grand total heat load QjOTAL * 

23 

^TCflAL ” ^ ^ ^TOTAL ^i 


The cooling load TR* expressed in tons is: 

TR " QtoTAL /1260'' (3.7) 

Assuming a safety factor of l.'l , the design tonnage 
TR is: 

TR = 1^1 Tr' (3.8) 


The cooling load can also be expressed as 

TR ~ c^ (h^ - h^) 3/12600 (3-9) 

where h. = specific enthalpy of air at exit (icj/kg- dry air), 

6 

hg = specific enthalpy of air at inlet (kJ/kg dry air) , 
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(b) Ambient and coach conditions; 

The calculation of ambient, coach inlet and coach exit 
condj-tions have already ueen discussed in sec# 2.'3# To 
determine the '.en^erature distribution inside the coach 
an iteration proeess is followed wherein at each section 
the following procedure is used: 

For any section <i*, we haves 


(a) For the 1st Iterations 


T^^ (3.10) 

At section 1 ( iie.,at inlet ) we haves 

^1 ~ ~ ^2 (S.'ll) 

Assuming the above, the values of are calculated and 
and used in the governing heat transfer equation f^r 
temperature distribution inside the coach to find the 
temperature T at section (i + 1)» Thus 


(b) For the succecdl-ig Iterations s 

A mean temperature, for section ’1* is 
calculated ass 


( 3 . 12 )' 




+ T_ )/2 

Using the abo-ve vailue of the new values of 
are calcu’ 

(i+l),‘ Hence, 


( 3 . 13)' 


U . are calculated to obtain the new temperature at section 
sj 


T 

^i+1 


T 


I 

i+1 


(3.14.) 
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The Iteration continues until the following constraint 
\+l - % < 0* 01 is s -tisfied* 


In this abo .e manner, the values of T„ are calculated 

^i 

for all the sections* calculate the humidity variation 
inside the coach, we use the governing equation of mass 
transfer. The rise in temperature ^ inside the 

coach is given hy 

A T = Tg - , (3.15) 

TEhe average temperature Tg_^ (‘^J) inside the coach is 
given by 


The relative humidity of air, 
given by 


$4 




where (p^)^ = 


P ^i 

0. 622 + w^ 


at any section *i* is | 

(3^7) 

(3.18) 


(c) Water requirement for the humidifier: 

Ullhe water requirement ( ) for the humidifier In 

kgA can he computed as follows: 

\ \ (wg - w^) (3.19) 


To calculate the air washer dimensions we use the following 
procedure; 

Assuming a face velocity of air through 

the air washer one gets: 



(3.21) 


whqr.e , _D ^ hvmldltler or jlr washer efficiency* 

and Z - JL-1 §J1- 

\ 

With h-0 ’5= mase transfer coefficient .(kgwA-Ja^-(kgwAg 

dry air))* 

“ volume of the air washer (m^) • 

ctW 

From Eq* (3.20) we haves 

Z - log^ ( 1 - n ) (3t22) 

From Eq,s. (3.21) and (3.22) we haves 

= - \ < 1 - l)/\ (3.23> 

For the outside conditions of air, the specific volume of 
moist air v(m^/kg dry air) is given by 

V = 287.2 (T,^ + 273.16)/[(P - Py ) 3c 10^] (3.'^) 

If is the volume flow rate of air (m^/h) through the air 
washer then 

(3.2^) 

To calculate the required length Cm) and face area 
(m^> of the air washer we use: 

Th^s the above quantities are calculated as: 


(3V26) 
13.. 27) 



raiHOUT KBGENEBATION: 


X. 


1 = 0; , 

T_ + 


T = 
aY 


with (p^)^ 


2 

PWj_ 






0.-622 + Wj^ 


WriH HEOEMlSAIIQIIi 


^1= Oi ^2 = ^:ci. ^1= 

(TJ„ + (li ), 


T =! 
•^aT 


e n 


*2 'n 




with [(PTr')i]n 


<4>n= 

I* ^"iV 
0,‘622 + 


Cd) SHCB AND USIKJRs 

Ihe performance of the humidifier at various places 
with different design conditions can be compared by use 
of the performance indices [ 2)i|.], * Sensible Bfeat 
Conversion Ifete ( SECR) * and *Unit Sensible Heat Conversion 
Rate (USBCR)t defined as; 


SHCR = fc^, C c^, 
and USHCR == SHCR/T^^^j^ 



) (kJ/h) 

(lcJy^-®0 


(3# 28) 
( 3 . 29 ) 


where o ” specific heat of moist air at inlet of 

pmi 

humidifier ( krAg - 




= spscific heat of moist air at exit of humidifier 

Ckj/kg- 

= temperature uf air at inlet of humiditler (°C)o 
= temperature of air at exit of humidifier (°C)* 


'^wbd ” "wet-bulb depression (°C )* 

The specific heat c (kJ/kg ~ ^C) of moist air is gi-ven by 

pm 

“pm = “pa » <3*3°^ 

For Cpg^ = 1, OOlf, ve have 

®pm = + l.'8%w (3.31) 

WITHOUT BSGEl^RATIONj 


“ 1 * 004 + 1 . 884 

♦ 


( 3 . 32 ) 

and c = 1,004 + 1. 884W5 

piu^ ^ 

• 

It 

C\J 

(3.33) 

WITH REGEHEIIA.TION; 

Vi = 3- 00“* + 884 

« 

> 

\ 

(3. 34) 

Cpjjj ~ 1,004 + l#8Sf W2 
^ 2 

• 

^b ~ ^2 

(3o3^) 

3,1) INPUT PiRAMETERS: 

The input data to the computer program for various • 

places on different routes is 

shown 

in Table 3«‘1» 




TABIiB 3.1 

LOKGITanES AND USTJt^JMS OF BIPCBTAin: STATIQUS (H 

DIFFERED ROUTES 


1 

P 

L 

A 

c 

E 

L 

T 

I 

T n i 
U(°) 

D 

E 

OUTSIDB 

DESIGN 

COI'IDITIO.. 

NS 

\ 

<°) 

(“) 

IiONGITUDE 
(EAST ) 

■ 

(*) 

(®C) I 

^wb 

(°c) 

SE Delhi 

29.0 

If 0, If 

23.9 

135.0 

45.0 

77,0 

12.0 

SE Kanpur 

27.0 

41.2 

21.5 

135.0 

45.0 

80,0 

30.0 

SS 

Ibllahahad 

25.0 

hX,7 

25.0 

135.0 

45.0 

81,0 

53.0 

SE Patna 

26.0 

37.9 

26.7 

135.0 

45.0 

85.0 

1.-0 

SE 

Jamshedpur 

23.0 

39.4 

28.4 

135.0 

45.0 

87.0 

1.0 

S 

Delhi 

29.0 

40.4 

23.9 

90.0 

90,0 

77.0 

12.0 

s 

Chandigarh 

31.0 

40.1 

23.9 

90,0 

90.0 

77.-0 

10^0 

s 

Bhopal 

23. 0 

40.2 

23.' 2 

90.0 

90.0 

78.-0 

1.0 

s 

Nagpur 

21.0 

42.6 

24.5 

90.0 

90.0 

79.0 

3.0 

).s 

Hyderabad 

17.0 

39.5 

26, 6 

90.0 

90.0 

77-.-0 

55.0 

L. SW Delhi 

29.0 

40.4 

23.9 

l^5.o 

135.0 

■“q 

12,0 

2i*SW Jodhpur 

26,0 

40.8 

25.3 

45.0 

135.0 

73.0 

1.0 

5. SW Jaipur 

27.0 

40.8 

21.7 

45.0 

135.0 

76.0 

30.0 

f. SW iliimedahad 

23.0 

41.-5 

29.0 

45.0 

135.0 

72.-0 

38.-0 

?, SW Baroda 

22.0 

40.3 

29.1 

45.0 

135.0 

72.0 
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The outside design condlt5.ons 
Table 3*1 are for the summer month of May,' The longitude of 
the standard meridiam is taken as 82 30* B, The 

declination angle and the equation of time 
have been obtained for May 21, The calculations for hour 
angle has been done for Indian Standard Time (X,S,T,) 
of 12,00 noon. The geometric dimensions of the coach have 
been taken from Appendix. The coach dimensions iirere 
obtained from drawing supplied by HDSO, Lucknow. Sius 
the final input values are: 

d = 20^, etjuatlon of time =* 3™ 37^» XST«12:00 noon 

“^^^ = 0.50, <«2 = 0.'40, •’gr = 

A = 1103 w/m®, B = 0-.190 C = 0.121 

p = 1,0132 bar, and rj = 0,75» 

The values of r.^ for glass window have been 
calculated as detailed in Appendix o* effect of 

regeneration on evaporative cooling has been studied 
by taking the following values of ® and r s 

® 3 0,-60, 0.70, 0,80 

r = 0 , 50 , 0,65, 0,80, 

The effect of moisture release has been studied 
by taking values of Mp as 0 , 050 , 0,100, 0,150 and 0,'200 
kg of molsture/kg dry air. 



3,2 COMPUIER PROGRAM: 


The computer program for the various calculations of 
both non-regnerative and regenerative cooling has been 
given in Appendix p • The function subprograms used are 
TIME, AG, PY, PS, Aar, FSS, FSG, REEL and the subroutines 
are REWON, FARGER and WETTEM, TIME returns the value of 
1 if L.A*T, is forenoon, 2 if IvA*T, is afternoon, 3 if 
UArT,' is noon. Depending on this value, the values of 
are calculated by the functions of program iG, PY calculates 
the partial pressure of water vapour in bars using Garriers 
equation and PS calculates the value of <T) in bars 
AT calculates the value of while FSS and FSG calculate 

the values of (F ) . and CF__)_j respectively. REFL 

S3 S J ' Sg S J 

returns the value of r^g for glass while NEWTON is a 
subroutine to calculate ( ^vi^si Newton -Raphson method 
to an accuracy of 1^. The subroutine PANGER gives the 
vaJ-ue of Convective heat '-ransfer coefficient over the human 
'':)dy from which comfort velocity V (m/s) is determined, 
rina.lly, the subroutine WETTEM uses Newton -Raphson iteration 
to find T^^ for the regenerative evaporative cooling system,' 



CmPTER 4 

ElStJLTS Aim DISCUSSION 

4*,1 EBTECT OF VARIOUS PARAMETERS: 

Effect of air flow rate on IR, AT and (without 
regeneration): 

Figures 4*1 and 4*'2 show the effect of air flow rate 
on TR and ^T, It is seen that with an incre4se in air flow 
rate, the tonnage increases by about Bfo due to a decrease 
in inside temperature of the coach as for higher mass 

flow rates is less.' With increase in the Inside heat 
transfer coefficient for free convection or forced convec- 
tion, the tonn^e TR and increase due to greater heat 
transfer to the coach. It is, however, found that free 
convection gives higher tonnage than forced convection due 
to small velocity. Hence the temperature Inside the coach is 
less for forced convection as cesapared to the free convection 
for the same inlet temperature of air* 

Figures 4*3, 4*4 and 4*5 exhibit the effects of air 
flow rate on the average teaperature inside the coach*' The 
results have been presented for various places of different 
routes at their respective design ambtait conditlotxs for 
summer.' Ihe average temperature Inside the coacbC ) 
decreases Int^riably with Increase in air flow rate* The 
ten^jeratures inside the coach are in the range -of 30 35^0 

which la not very far from the comfort conditions. 
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4*Lr2 Iffeot of train yelocity on TRr 

I^gure 4«6 shows ihe effect of train velocity on TR, 

It decreases by about 1 to 3 ^ with an increase in train 
speed from 60 to 120 IxmAi* This is due to the fact that 
with higher train velocity, h^ increases rapidly than the 
overall heat -transfer coefficient U due to which the sol-air 
temperature gets decreased. Hence their net effect causes 
decrease in TR, The inside temperature at higher velocity 
is found to be lower than that at lower velocity of the 
train# This has been observed for the regenerative case 
as well, 

4.1«‘3 Effect of regeneration: 

a) Effect of recirculation: 

Figures 4*9 to 4rll show the effect of recirculation 
on TR, A T and Tg^^‘ Increased recirculation renders the 
lower tonperature inside the coach as compared to the non- 
regenerative case* The tonnage for regenerative cooling 
is therefore higher than that of the non -re generative 
case. Further, the higher r:; circulation causes better 
comfort conditions inside the coach#' 

b) Effect of heat exchanger effectiveness: 

The higher heat -exchanger effectiveness also 
produces the same effect as increased recirculation as shown 
In Pigs# 4*12 to 4.14. Here two temperatures inside 
the coach are lower than that of the non -regenerative case# 
However, It is found that increasing recirculation is more 
effective in bringing down the temperatures inside the coach 
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than increasing heat exchanger effectiveness. This is 
supported by the f^Gt ^hat recirculation decreases the 

i 

temperature ly about '4»7^ ■where s heat exchanger effective- 
ness, decreases the same by only 3%* 

4«1*4 Effect of moisture release: 

, Jigures 4* 15 and 4*16 show the effect of moisture 
'release iori tonnage and temperature distribution inside the 
coach. The higher the moisture release, the lower the 
temperature inside the coacho This is due to the fact 
that the sensible heat of air is partly used by moisture 
for its latent heat of vapoii.zation causing decrease in 
the air temperature. Bat the tonnnage increases by aTbout 
2* 3% as the coach temperature is lowered. 

4* 1*5 Effect of air flow rate on ft : 

1d.th increased air flow, rate the water requirement 
for the humidifier increases linearly as shown in Tigs* 

4.17 and 4* 18* toiuporaturo at inlet to humidifier 

decreases by about 0.'2°C depending upon the rate of 
recirculation. For the same evaporative cooling efficiency 
(n = 0,75 ), the specif: c humidity at exit of humidifier 
is less compared to that without regeneration. Hence ,the 
value of is less for the regenerative process.' It is 

to be noted that the rejeneratlve process requires less' 
water supply and renders better inside conditions as 
compared to non -re generative case. 

4.1*6 Effect of air flow rate on SIKE and USIKRt 

Figure 4«il9 shows that the performance indices SHOE 
and USHCE increase linearly with air flow rate. The water 



requirement tlie Yalues of SHCR and USHCR for various 

places with given desigr ambient conditions for sunaner are 
presented in Table 4.1 for an air flow rate of 10,000 kg/h. 

It is observed from Table 4. 1 that the indices of performance 
SH5R and USHCR decrease with increase in the relative 
humidity* df the ambient air. So , for effective evaporative 
cooling, the relative humidity of the ambient air should be 
as low as possible. To get comfort conditions inside the 
coach the air circulation rate should not be less than 
10,000 kg An, However, increase in air flow rate would 
Increase the space required for the humidifier, the size of 
the blower and other components.' Therefore, it is suggested 
that the air supply should be kept in the range of 10,000 to 
1?,000 kg/h. 

The computed values of comfort conditions inside the 
coach for various mas £3 fio# rates with regeneration 
( e = 0.60, r = 0 , 50 ) for Kanpur are presented in Table 4.2.' 

The required a*'* r-washer dimensiois are also Indicated for 
a face velocity of 2.4 m/s and quantity hjj taken to be 
8000 kgwAi-m^- (kgw/kg dry air). It is easily seen from 
Tables 4* I ahd 4. 2, that regeneration gives lower inlet and 
exit temperature of the coach compared to the non -regenerative 
case for Kanpur." The humidifier dimensions lie within the 
Space being available underneath a coach ( as per the RDS0*s 
infomationX 



TABLE lf .1 

VALUES OF aAND PSRFOBMANCE INDICES FOR VARIOUS PLACES 
(WITHOUT REGMERATION ) 


SIS 

|T{ ir 
B*® o 

E <8 
tt>+s o 
Eh fij O 


Eo ^ 

^ O 
+> 4 > 

- 

Si n ^ 

&**rl O 
etJ 

o 
o 


m 

<a> 4 J 
Eh (\J 


oS 

..C 


B 


. h> 

{D v-*' 


«■§ 

B !3 




c 

o 

o 

p 

to 

*H 

xn 


a 


*G 4 


•g 0 

O 


•P 

in c 

(D O 

*H ^ 

,0 Q 

■§:s^”’§ 


0 ) 

H 


<D 

-P 

El 

O 


£ 

ca 


H r-i Vv \jf\ H ^ j:j' VO ^ Jt 

^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

^^<^ro(vi<^rororo<vioorornro 


O^OJT^CvlO'ii^.OCO CMxrxH O 

9R'5pi5>^'^c)Ocva^cNCK\DC^f-l 
CVJCVlOJOJmCMCM CJCU 


rr>oo ^3 

'tTN^VO VO 

-±lf\vOi^mvo Jr \ts ^ ^ ^ 

^ ^ #* •v ^ <1^ ^ 


H H H 


H H H 


- ha 1 

-dh ^ 

•» <jv CO 

H HI 


?N ^ jdh V\ H ^H NCOOO cry CT^ no 

H H CM \r\ Xf\ O CO vO O OC? CN CN \f\ 

Its \0 \r\ cr\ 00 \r\ \f\ Tjpv jj* j- \f\ (r^ fry 


^£\ ^ \p CVI^vDmcMN 
OJHCMJ-5c\JCM<MroCMH-ih 

^ «k 

oooooooooo 


gfS 3 - J: 


000 


OV ITN O 

m H \r\ M 3 00 
CM CM CM CM CM 


lN^ON<M\r\vO rniNO H 


CM cJ ^ 


lf\ 

CM 


■\ ^ cA 

M CM CM 


On 

CM 


y^cMiN^ON^H cMvo\r\cOcx)xrscn 


O H H 
J* -dh -d" 


05 

^ nJ 




^ •? 
S i 

-d «Mj 

03 ♦H 



^ lu w -n r— 1 f-i 

'm S'ji;? 9 3 5 ^-p-D'H 

iSigg«>«„B,gg 


‘ ? /-Jf ^ ' f ^ ♦ ' ♦ O H 

rrlCMfOJ-VvvOiv.CX)OSr4H 


B 


k9 


50 


TABia; 4*2 

COJiPH CONDITIONS AND AIR«WASHEE DIMENSIONS FOR KANFUR 
(Tdb = °C, = 21,5 °C) WITH BEGENEB4EI0N 


ft 

'se. 

— 

Vb 

5 

a 

r) 

~ — r 


’f 

^2 

T 

e 

$ 

7 

cv 

^aw 

■^aw 

.®c) 

L°0) 



(kg/h) 

' (®C) 


(®c) 

(m/s) (m) 

(m^) 

i 8.8 

20.8 

0 V 18 

0.75 

5000 

25.3 

0.75 

33.3 

0.50 

0, 88 

1,’67 

0*52 

(7.7 

2D.-5 

0.20 

0.75 

10,000 

24.' 8 

0.75 

29.1 f 

0,60 

0.78 

1.-68 

1,03 

17.2 

20.3 

0.20 

0.75 

15,000 

3 f . 5 

0.77 

27.9 

0,-65 

0.-35 

1,68 

1.55 

l7,0 

20,2 

0.-20 

0.75 

20,000 

af.4 

0,78 

27.0 

0.68 

0 V 3» 

1.69 

2,06 


Th® water -reQulTQiiidnt for the humidifier for Kanpur assuining 
i servio© duration of 6 hours is found, to be 370 litres with 
30 regeneration, whereas for a regenerative evaporative cooling 
astern (e « 0,60, r = 0,'50),the same, comes out to be only 305 
Litres of water. Hence the water consumption is less for regenerative 
system in addition to the advantage of better comfort conditions 
Lnslde the ooach. 


I 
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W2 IXPERIMSNTAL EESULTS: 

To study the temperature distribution Inside the coach 
experimentally, a coach model was fabricated from aluminium 
sheet with a ratio of 1:20 compared to the dimensions of 
the prototype railway coach* The model was kept inside a wind 
tunnel and hot air was blown across it. The temperature 
of the hot air was maintained at the actual ambient con- 
ditions of Kanpur during summer ', = ^l^'S^C), The air 

velocity inside the wind tunnel was used to simulate the 
actual moving train. The air flow rate through the coach 
was measured by calibrated rotameter, A heating coil was 
used to sinmlate the heat release from persons. The 
experimental set-up is shown in Fig, 4,20, Four calibrated 
thermocouples were used to measure the temperatures of air 
at various points inside the coach and the ambient air 
temperature. The thermocouples were connected to the direct 
temperature indicator potentiometer. The results are ^own 
in Fig, 4*21 where the dimensions of the experimental set.hp 
is also shown. The variation in temperature Inside the coach 
is linear which is similar to the predicted nature of 
temperature variation.' However, the tempersiture Increases 
sharpl 3 ^ in the initial section of the coach. This Is 
due to the fact that the temperature of the supply air at 
inlet to the coach is about 22*^0 which increases sharply 


duel, to the heat release mainly because of high thenKil 
eooauetivlty t.f the aluminium coach mod^, 


rt) 

h ma mm sm WM 


I 
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If, 3 PROPOSED SXSM: 

Figure if, 22 gives as schematic representation of the 
proposed evaporative cooling system for the actual rail^^ay 
coach. Water for the evaporative cooling system is taken 
from the main supply line at the station, The stored water 
is pumped through the spray nozzles of the humidifier and the 
humidified cooled air is passed through the eliminator 
to remove dust and water droplets,' This air is then 
introduced into the coach to produce comfort conditions,' 

In the actual case, to maintain better comfort conditions, 
the humidified air should be supplied to passengers 
directly through grilles. This analysis may be regarded 
as the extension of the p3?esent work. 
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O.) 



Fig.4.3. Variation in Tov with air flow rata tor different places 
without regeneration. 



H 



Fig.<.4.Voriation in Tq^ with air flow rote for ditf«r«nl phico* 
without r»g«o«ration. 



TRAIN VELOCITY = 60 km/h 



itIq (kg/h) 


Fig.4.5. Voriation in Tqv with air flow rate tor different places 

without regeneration. 



place: DELHI 

route; south-east 
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Fig.4.6.Voriotion in tonnage with train velocity. 
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Fig, 4.11. Temperature distribution iriside the coach for fa) without regeneration (h) with re- 
generation (Heatexchanger effectiveness (€) =0.60) (Place- Kanpur). 



rrt« convection at sides and floor 
Forced convection at oil sides {Vfc*aS m/s > 
train velocity a 60 km/h 

i £ » 0.80 '1 
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Fig. 4 12. VariQjion in tonnage with air flow rate for ia) without 
regeneration end ib) with regeneration using 50% re- 
cifcuiution of exit air from the cooch {Place -Kanpur). 
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4.17. Variation in with air flow rat® for ia) without regenera- 
tion and (b) with regeneration {Heat -exchanger cffcctivencsi 
* 0.60) (Place • Kanpur) 
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Variation in with air flow rate for to) without 
ntrotion and fb) with regeneration using 50% recircu- 
lation of exit air from the coach CPloce - Kanpur ) 
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Rg. 4.!a Variatii>n in SHCR and USHCR with air flow rate for 
(q) without rogancration (b) with rcgcncrotion (CsO.60. 
r»a50) (PUCE -KANPUR). 
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Fig. 4.22. Schematic representation of proposed evaporative cooling system 




5.1 CONCUJSIOHS: 

yi'om the present study the foUowlng conclusions are 

arrived at t 

1* ^he cQjnputer program has been developed to calculate 
cooling load, temperature distribution Inside an 
evaporatlvelj -cooled railwav coach and for non- 
regeneratlve and regenerative evaporative cooling 
systems.' 

2. The results have been obtained for a wide range of 
operating variables, for three different train routes 
(SE, S and SW ) , 

3# Tho cooling load TR decreases by 1 to 3?^ with an 
increase in train velocity from 60 to 120 toaAv 

4* Higher air flow rates give better comfort conditions 
within the coach But the space req.ulrement for 
the humidifier would be large in addition to 

bigger size of the blower and other components. 

Hence, an air flow rate in the range of 10,000 to 
15,000 kg A is recommended in view of the availability 
of for the same underneath the coa^h. 

5v The regenerative evaporative coollpg systi^m Is 

advantage's Compared to non -regenerative evaporative 
dopllng iti terms of achievement of better 

ik»ido ©laidiiions less water consumption, 
fbr a six hourly sqplce of th^ system per day, the 

i 
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regenerative evaporative cooling system requires 
only 305 litres of water compared to non -regenerative 
system where the requirement is 370 litres of water*' 

6* !I!he performance of the regenerative cooling system 
can be enhanced by either increasing the recirculai . - 
on rate of exit air on the he at -exchanger effective- 
ness or both. However, it is found that the former 
is more effective than the latter, because the 
temperature inside the coach decreases by about 
4*7?^ in the former as compared to only in the 
latter. 

7. The computed dimensions of the humidifier can 
easily be accommodated underneath the coach, thus 
making evaporative cooling practically feasible for 
a prototsrpe railway coach. 

8. The tonnage of the evaporative cooling system 
comes out to be in the range of tons which 

is only half of that of the conventional 10 ton mecha- 
nically refrigerated air conditioning system 
used for railway coaches.' 

5. *2 SCOPE FOa FIMJKE WORK 

There is considerable scope for the exteniiion of 
the present work. They are: 

1,' Analysis of heat and mails H-ansfer folp cool;)Lng 
. load calculations of the coach can be Carrledl out 
when the air is'' Introduced through grilles fa*om 
both ends of the coach;' 
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2# The effect of activity of occupants inside the coach 

' , 1 ’ll 

on the cooling load can be. studied* 

3* !IhG cooling load may be calculated using combined 
free and forced convection for inside heat transfer 
coefficient. 


4. The sizes of the humidifier, regeneraUve heat 

,, ' , .ii '■*'-' ' • ' 

oxohangor, supply duct, blowers and pumps can 
be optimised on the basis of overall cost and 


energy requirement; 

5; An experimental investigation can be undertaken 
for a prototype railway coach.’ 
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BVALUATIOU OF COACH PARAMETERS 


A.1 EVALUATION OF COACH PERIKETER (p ) 

c 

We have the dimensions of the coach as Lj^ = 21,'337 m, 

= 3 * 06 l m, = 2,075 m,- 

= 2I^ + V CA1> 

where = perimeter of the curved ceiling of the coach (m) 

Pg perimeter of the other three sides of the coach (m) . 

To calculate P^ we consider the ceiling as approximately 
elliptical with major axis as ’a* and minor axis as 
Then equation of the ellipse is 



Putting X = ae. sin 0 , y =.b cos 0 

dx = a cos 0, dy = - b sin 0 (A3) 


?! = Jytdx)^ + (dy) 

•A „ 




cos^0 + b^ sin^0 d0 - 
a /jTrTli?7a?) sin^0 d© 


Putting = IJy/ar f in Eq, ( A 1+ ) one gets! 

jjy^ n il ' 

P. = a 1 - k^ sin^e ) d0 
'^o 

using the following Eq, (A 6) from [ 2?] ? 


(A if ) 


(A 5) 




\/l - k^sin^0 d© -(2/tc) + sin § cos § [l~ , - k^ + 

^ k^A. 


^ A + 

2. if ^ 

E($,k) 


2, if. 6 
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7t/2 






where 

E = 

r j 

V/ V 

0 

' 1 - 

sin^e de 

= l-[i - 

. (i_ 
^2 





- ( 

2.1^6 ^ 

j£ 

5 

...] 

if k^ <1 

For § 

= Tt we have 







}1 

V i 

1 - 

sin^e . 

der = 2 B 

orPj_ = 

2aS5 

(A 7) 

(A. 2) 

krei. 

of cross-section (A_) 





A 

*c 

it 

h 




(A 8) 


2 

= area of the semi -ellipse under the ceiling (m ) 
Ag s= area bounded by the walls and floor in?) 


For send. -ellipse 



s : 3 , t )/2 


(A 9> 

and 



(A 10) 


The final results 

obtained are: 



p.^ = 10,665 m 
^ = 7.533 



appeijdix b 


BVA2.0AIIM OF T^; F« BEOEHBBATIOK 


From: Eq, ( 2 , 83 ) 


P (t'v.) ** 
S wb^ 


e 


M 


(B 1) 




M = [ 7.2137 + 1.1520 X 10-5 X 10-9 

<i;4l3 + 273.16) .^10)2] (.. Ay tjl 1 1](B 2ii 

1^^+273.16 

. , [ 1.0132 -P,(J;,;>][ la- 

CPv^ “ l>s ^ - T"^ 3) 


15W - l-% T, 


vb 


f.™ S - , (T [^•'>132 - 1 ^ . T 

~ “ 


ualng Nbwton-Baphson metv.^d 


Vl 


^^wb^ k 


151*7 - 1.1(4 


I ,_ I 


.p_ 0 




(B If) 


with k = no* of iterations 

f 

t 

wb 


where f* (T„J) « 


( 151+7 J..-l|4 ) ( Ta - ’^wb H221.2 e-l* _dil),/(151+7 J..l(l+(T„t^ ) 

dT. 


’wb 


_ ^^a“^wb )CXfp)'3(2.^2i* 2 e )/(1.5lf7-d.*4^(T,^^ ) 

Putting (T^)i*= T^Tjj above iteration is performed k tlJiies.' 
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jgmpix cr 

BmUAIION OF H^LEOTAIilCB OF GLASS 


A subroutine for calculating the reflectance of glass 
been incorporated in the main program which uses 
Lagran^ian interpolation method for calculating the value of 

®CS 


Ln(x) 


Ij^ (x) 


where x^,.. ..x^ 


not necessarily equally spaced and 


1 q(x) = (x-Xj^) (x-Xg) ......... (x-Xj^) 

Ij^(x) = (x^Xq) Cx-Xj^^j^). ... (x-Xj^) ,h=l, 

• • * n»JL 

lj^<x) = C 3 |;-Xq) (x-Xj_) ... .. 

The data used for interpolation taken from graph [ 3 j ] is 
given as under: 


TABLE B*1 

VALTJES OF BEFLECTAECE FOR VARIOUS INCIDIMT AHjLES FOR GLASS 



r ^ 

cs 

50 

0.06 1 

60 

0.09 

70 

0. 165 

80 

0,38 i 

i' 

90 

1.00 1 

f 


%if 

®cs linear interpolation is used: 

r^g = 0,-02 + ( 8 X 10**^) e^g (C 1) 

Bie vindow glass is of double strength quality with the 
following values: 

tg = 3.175 
K = 7.64 I 10“2 
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COMPUm PROGRAM 
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